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EFFECT  OF  UNCONFINED  LOADING 
ON  THE  UNFROZEN  WATER  CONTENT 
OF  MANCHESTER  SILT 


J.L.  Oliphant,  A.R.  Tice  and  R.  Berg 


INTRODUCTION 

One  variable  that  is  basic  to  understanding  the 
frost  susceptibility  and  strength  characteristics  of  a 
partially  frozen  soil  is  the  unfrozen  water  content 
Wu.  The  unfrozen  water  content  has  been  measured 
directly  using  calorimetry  (Williams  1963,  Ander¬ 
son  1966)  and  nuclear  magnetic  resonance  (NMR) 
(Tice  et  al.  1978).  Thermodynamic  considerations 
have  made  it  possible  to  estimate  IV u  from  soil  mois¬ 
ture  characteristic  curves  (Williams  1964,  Koopmans 
and  Miller  1966)  and  from  other  data  relating  the 
activity  of  the  soil  water  to  the  water  content  (Low 
et  al.  1968). 

The  effect  of  overburden  pressure,  or  surcharge 
pressure,  on  the  unfrozen  water  content  has  also 
been  predicted  using  thermodynamic  arguments 
(Edlefson  and  Anderson  1943,  Williams  1964,  Low 
et  al.  1968).  However,  direct  measurements  of  the 
effect  of  pressure  on  unfrozen  water  content  are 
few.  A  small  amount  of  data  for  which  the  pressure 
effect  is  significant  has  been  reported  (Tsytovitch 
1959).  Other  experiments  that  show  migration  of 
ice  in  the  soil  from  areas  of  high  pressure  to  areas 
of  low  pressure  (Vialov  1959)  also  demonstrate  a 
significant  pressure  effect. 

In  this  study  the  unfrozen  water  content  of  a 
silt  from  Manchester,  New  Hampshire,  was  deter¬ 
mined  as  a  function  of  total  water  content,  temper¬ 
ature  and  surcharge  pressure  using  the  NMR  tech¬ 
nique.  The  data  were  analyzed  statistically  to  de¬ 
termine  relationships  between  the  unfrozen  water 
content  and  each  of  these  variables.  The  effect  of 


surcharge  pressure  was  then  compared  with  that  pre¬ 
dicted  thermodynamically  by  various  forms  of  the 
Clausius-Clapeyron  equation. 

MATERIALS  AND  METHODS 

The  soil  used  in  this  investigation  was  sampled 
from  the  vicinity  of  Manchester,  N.H.,  and  is  desig¬ 
nated  Manchester  or  New  Hampshire  silt.  The  phys¬ 
ical  characteristics  of  the  soil  have  been  well  defined, 
and  numerous  tests  of  its  freezing  behavior  have  been 
performed  (Koopmans  and  Miller  1960,  McGaw  1972, 
1973).  Six  200-g  aliquots  of  dry  soil  were  mixed 
with  distilled  water  to  form  mixtures  containing 
0.10-0.225  gof  H20  per  g  of  soil.  The  soil-water 
mixtures  were  sealed  and  allowed  to  equilibrate  for 
one  week.  Following  moisture  equilibration,  eight 
samples  at  each  water  content  were  compacted  in 
three  layers  each  in  a  1.57-cm-diameter  by  4.02-cm- 
high  mold.  The  compaction  hammer  weighing  300.9 
g  had  a  free  fall  of  1 1 .91  cm  and  exerted  a  compac- 
tive  effort  of  3.5  x  106  ergs  per  blow.  Seventy-five 
blows  were  given  to  each  sample.  Following  com¬ 
paction,  each  sample  was  trimmed  to  3.2  cm  in 
length,  and  a  hole  0.5  mm  in  diameter  was  drilled 
along  its  axis  to  accommodate  a  copper-constantan 
thermocouple.  Each  sample  was  inserted  in  a  1 .68- 
cm-i.d.  test  tube  and  sealed  with  a  rubber  stopper 
to  prevent  moisture  loss.  The  test  tubes  containing 
the  compacted  samples  were  placed  in  a  temperature- 
regulated  bath  containing  a  mixture  of  ethylene 
glycol  and  water  at  a  temperature  of  about  +0.75°C. 


The  coolant  in  the  bath  was  vigorously  circulated 
with  two  submersible  pumps,  and  the  temperature 
was  controlled  with  a  Bayley  proportional  tempera¬ 
ture  controller. 

A  Praxis  model  PR-103  pulsed  NMR  analyzer  was 
operated  in  the  90°  mode  with  a  0.1  -s  clock  and  at 
a  fast-scan  speed.  The  first  pulse  amplitude  in  the 
90°  mode  was  measured  for  each  sample  starting  at 
about  +0.75°C.  The  test  tubes  were  sequentially 
removed  front  the  bath,  wiped  dry  and  inserted  in 
the  NMR  anrlyzer.  After  about  4  s  (the  time  re¬ 
quired  to  record  the  sample  temperature  and  the 
NMR  pulse  amplitude),  the  samples  were  reinserted 
in  the  bath.  When  all  the  samples  were  analyzed, 
the  temperature  was  lowered  about  0.2°C  and  the 
measurements  were  repeated.  For  each  sample  four 
readings  above  freezing  were  used  to  form  a  ratio 
between  the  NMR  signal  intensity  and  the  sample 
water  content. 

Following  the  last  observation  at  above-freezing 
temperatures,  the  test  tubes  containing  the  samples 
were  placed  in  a  coldroom  at  -20°C  and  allowed  to 
freeze  overnight.  The  next  morning,  circular  lead 
weights,  which  were  also  cooled  to  -20°C,  were  in¬ 
serted  in  the  test  tubes  and  brought  to  rest  on  two 
1 ,25-cm-thick  grooved  Teflon  disks  that  were  placed 
on  the  soil  (Fig.  1).  The  grooved  disks  ensured  that 
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Figure  1.  Schematic  representa¬ 
tion  of  a  test  tube  containing  a 
frozen  soil  plug  with  a  lead  weight 
surcharge. 


the  thermocouples  did  not  interfere  with  the  loading 
process  and  acted  as  thermal  insulation  between  the 
soil  and  the  weights.  Six  lead  weights  ranging  front 
75  g  to  I  POO  g  were  added  to  six  samples  at  each 
water  content.  Two  samples  at  each  water  content 
were  unweighted  to  provide  a  check  on  the  reproduc¬ 
ibility  of  the  measurements.  The  unweighted  samples 
and  the  two  sets  of  samples  with  smaller  w-eights  were 
sealed  with  rubber  stoppers  to  prevent  moisture  losses. 
The  larger  weights  extended  above  the  tops  of  the 
test  tubes  and  were  sealed  with  a  thin,  circular  finger 
cot  with  the  end  cut  off  (Fig.  1 ). 

The  samples  were  removed  from  the  coldroom 
and  quickly  placed  in  the  temperature  bath,  preset 
to  about  -13.5°C.  The  top  centimeter  of  the  test 
tube  and  the  larger  weights  extended  above  the  top 
of  the  bath  (Fig.  2),  which  caused  two  difficulties. 

First,  one  set  of  weighted  samples  was  ruined  be¬ 
cause  moisture  accumulated  on  the  cold  weights  and 
then  permeated  into  the  soil.  This  problem  was  solved 
by  installing  a  dehumidifier  in  the  room  and  by  blow¬ 
ing  air  over  the  samples  with  a  large  fan. 

Second,  the  lead  surcharge  weights  that  extended 
above  the  bath  cover  conducted  heat  into  the  soil 
samples.  This  became  apparent  when  the  recorded 
temperatures  for  the  samples  with  the  heaviest  weights 
exceeded  0°C  and  when  the  NMR  data  did  not  ap¬ 
proach  the  value  expected  for  a  thawed  sample.  To 
determine  the  cause  of  this  discrepancy ,  which  was 
thought  to  be  related  to  heat  conduction  through  the 
lead  weights,  the  following  test  was  conducted.  The 
temperature  of  the  bath  surrounding  the  samples  was 
set  at  about  0°C,  and  the  samples  were  allowed  to 
reach  thermal  equilibrium.  Temperature  readings 
were  taken  on  all  samples.  The  weights  were  removed, 
and  the  depth  of  each  thermocouple  below  the  soil 
surface  was  measured.  Depths  ranged  from  0.5  to 
1 .6  cm  and  averaged  1 .04  cm.  All  the  thermocouple 
holes  were  then  redrilled  to  a  depth  of  2.7  cm.  and 
the  thermocouples  were  reinstalled.  The  2.7-cm-depth 
corresponds  to  the  limit  of  sensitivity  of  the  NMR 
equipment.  The  weights  were  then  reinserted,  and 
the  sample  and  weights  were  again  allowed  to  reach 
thermal  equilibrium.  Temperatures  were  again  re¬ 
corded  for  all  the  samples,  and  this  time  the  temper¬ 
atures  were  within  ±0.003°C.  indicating  that  the 
temperature  perturbation  did  not  extend  into  the 
region  of  interest,  that  is,  the  region  surrounded  by 
the  NMR  detector  coil.  The  temperatures  were  also 
within  ±0.003°C  of  those  measured  near  the  top  of 
the  samples  when  the  small  weights  were  used.  There¬ 
fore,  when  the  data  were  analyzed,  the  temperatures 
of  the  samples  containing  the  heaviest  weights  were 
changed  to  those  recorded  for  unweighted  samples 
immediately  adjacent  to  them. 
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Figure  2.  Weighted  samples  in  the  constant  temperature  bath. 


Following  temperature  equilibration  at  -13.5°C, 
each  sample,  including  its  surcharge  weight,  was  se¬ 
quentially  removed  from  the  bath,  wiped  dry  and 
inserted  in  the  NMR  detector.  After  about  4  s  the 
samples  were  reinserted  in  the  bath.  When  all  the 
samples  were  analyzed,  the  bath  temperature  was 
changed  and  allowed  to  equilibrate.  The  process  was 
repeated  until  data  for  complete  phase  composition 
curves  were  obtained.  NMR  signal  intensities  were 
measured  on  each  sample  at  18  temperatures  between 
0°  and  -13.5°C.  There  were  eight  samples  at  each 
water  content  (six  weighted  and  two  unweighted), 
and  measurements  were  made  at  six  water  contents. 
This  resulted  in  144  readings  at  each  water  content 
and  864  for  the  entire  experiment.  After  the  final 
observations,  water  contents  were  determined  grav- 
imetrically,  and  a  ratio  of  sample  water  content  to 
thawed  first  pulse  amplitude  was  developed.  Unfro- 
water  contents  were  calculated  by  multiplying  first- 
pulse  amplitudes  by  the  respective  ratios  to  obtain 
values  of  the  unfrozen  water  content  at  each  tem¬ 


perature  (Tice  et  al.  1978).  The  unfrozen  water  con¬ 
tent  vs  temperature  data  are  listed  in  the  Appendix 
for  each  water  content  and  surcharge  weight. 


DATA  ANALYSIS 

According  to  the  Clausius-Clapeyron  equation 
(Edlefson  and  Anderson  1943,  Williams  1964)  if  the 
surcharge  on  a  frozen  soil  causes  the  pressure  on  the 
ice  and  water  phases  to  increase  by  the  same  amount, 
then  each  kilogram-force  per  square  centimeter  of 
pressure  will  shift  the  unfrozen  water  content  vs 
temperature  curve  by  0.0073  Celsius  degrees.  If  the 
surcharge  increases  only  the  ice  pressure  while  the 
water  pressure  remains  constant,  then  the  temperature 
drops  by  0.088  Celsius  degrees  for  each  kilogram- 
force  per  square  centimeter.  In  either  case  the  effect 
of  a  surcharge  is  to  shift  the  unfrozen  water  content 
vs  temperature  curve.  In  the  following  analysis,  then, 


ft 


any  surcharge  effect  is  assumed  to  have  an  equiva¬ 
lent  temperature  effect. 

The  surcharge  effect  on  unfrozen  water  content 
was  isolated  by  first  finding  an  empirical  equation 
that  would  reasonably  describe  the  unfrozen  water 
content  vs  temperature  curve.  The  parameters  in 
this  equation  were  determined  by  using  a  nonlinear 
least-squares  analysis  of  all  the  data  from  samples 
having  the  same  total  water  content.  In  this  initial 
analysis  the  differences  in  surcharge  among  the 
various  samples  were  ignored.  Then,  wherever  the 
independent  variable  temperature  T  occurred  in  the 
empirical  equation,  the  combination  of  independent 
variables  T and  surcharge  S  were  substituted  in  the 
form  T  +  KS,  where  K  is  the  constant  that  defines 
the  temperature  equivalence  of  the  surcharge  and 
has  units  of  °C-cm2/kgf.  (The  values  -0.0073  and 
-0.008°C-cm2/kgf  from  the  Clausius-Clapeyron 
equation  are  K  values.)  The  revised  empirical  equa¬ 
tion  was  then  fit  to  the  data,  again  using  the  non¬ 
linear  least-squares  analysis,  this  time  including  the 
effect  of  surcharge.  In  this  way  the  experimental 
value  for  K  could  be  determined.  Also,  it  could  be 
determined  if  including  surcharge  effects  significantly 
improved  the  fit  of  the  equation  to  the  data,  that  is, 
if  the  effects  of  surcharge  were  significant. 

An  empirical  equation  used  by  Anderson  and 
Tice  (1972)  to  relate  temperature  and  unfrozen  water 
content  is 

Wu  =  <x6»  (1) 


where 

lfu  =  unfrozen  water  content  (in  7<  of  dry 
soil  weight) 

0  =  temperature  below  freezing  (°C) 
a  and  (3  =  empirically  determined  parameters. 

Although  this  equation  fits  the  general  shape  of  It' 
vs  6  curves  quite  well,  the  poorest  fit  was  at  the  high¬ 
est  unfrozen  water  contents  (at  temperatures  slightly 
below  0°C)  on  the  Manchester  silt  samples.  Figures 
3  and  4  illustrate  this  for  a  typical  sample  (20%  total 
water  content  and  no  surcharge).  In  Figure  3  the 
data  and  the  best-fit  curve  are  shown.  In  Figure  4 
the  residuals  (the  differences  between  predicted  and 
measured  values  of  Wu)  are  shown  as  a  function  of 
temperature.  The  curve  for  eq  1  fits  especially  poorly 
in  the  region  between  0°C  and  about  -1°C. 

The  region  from  0°C  to  -1  °C  is  where  the  unfro¬ 
zen  water  content  changes  the  fastest.  Therefore, 
this  is  the  region  where  the  unfrozen  water  content 
is  most  sensitive  to  any  surcharge.  The  empirical 
equation 

Wu=A+B[l-exp(CO))IO  (2) 

fits  the  data  better  than  eq  1 ,  especially  at  tempera¬ 
tures  just  below  0°C.  In  this  equation.  A .  fland  C 
are  empirical  constants.  Equation  2  has  the  property 
that  at  0=0,  WU=A  +  B  (-T) 

The  parameters  in  eq  2  were  determined  at  each 
total  water  content  by  using  all  the  data  obtained 


Figure  3.  Unfrozen  water  content  of  a  Manchester  silt  containing  20%  of 
dry  weight  total  water  and  with  no  surcharge. 
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Figure  4.  Residuals  showing  goodness  of  fit  of  eq  1  and  eq  2  for  the  data 
shown  in  Figure  3. 


from  samples  having  that  water  content,  regardless 
of  surcharge.  Then  each  set  of  data  was  used  to  find 
the  parameters  in  the  modified  equation 

WU=A+B  |l  -exp  [a0+A'S)]}/(0+A'S) 

(3) 

taking  into  account  the  surcharge  5.  The  values  of 
A,B.C and  K  are  listed  in  Table  1  for  each  set  of 
samples.  Also  listed  are  the  values  for  Wu  when 
6  +  KS  =  0.  These  values  show  no  significant  increase 
with  increasing  total  water  content,  indicating  that 
the  amount  of  unfrozen  water  does  not  depend  on 
the  total  amount  of  ice  in  the  samples.  The  listed 
values  of  K  show  no  discernible  changes  as  a  function 
of  total  water  content.  The  six  values  have  a  mean 
of  -0.297°C-cm2/kgf  and  a  standard  deviation  of 
0.048°C-cm2/kgf.  The  values  for  K  that  were  found 
experimentally  are  considerably  higher  than  those 
predicted  using  the  Clausius-Clapeyron  equation,  i.e. 
-0.088  or  -0.0073°C-cm2/kgf. 

There  are  several  ways  to  test  the  hypothesis  that 
eq  3  fits  the  data  significantly  better  than  eq  2.  This 


is  equivalent  to  saying  that  surcharge  has  a  significant 
effect  on  unfrozen  water  content.  Because  the  values 
for  K  in  Table  1  are  relatively  consistent  and  because 
zero  is  not  within  two  standard  deviations  oi  the 
mean  value  of  A,  it  is  likely  that  the  surcharge  has  a 
significant  effect. 

The  sum  of  the  squares  of  the  residuals  obtained 
when  using  eqs  2  or  3  as  the  model  can  be  divided 
by  the  degrees  of  freedom  attributable  to  the  resid¬ 
uals  (141  for  eq  2  and  140  for  eq  3)  to  obtain  values 
of  the  residual  mean  square  (Table  2). 

A  ratio  or  he  residual  mean  square  values  can  be 
formed  to  estimate  the  F statistic.  At  the  90%  level 
of  confidence,  all  six  values  for  F  show  that  eq  3  fits 
the  data  significantly  better  than  eq  2  (Table  2).  The 
F statistic  is  rather  sensitive  to  the  assumption  that 
the  populations  being  considered  (in  this  case,  the 
residuals)  are  normally  distributed.  Even  though  eq 
2  fits  the  data  better  than  eq  1 ,  it  is  still  not  a  “cor¬ 
rect”  model  in  that  the  residuals,  as  shown  in  Figure 
4,  are  neither  randomly  nor  normally  distributed. 
Therefore,  using  the  F  statistic  to  compare  dispersion 
is  probably  not  valid. 


Table  1 .  Parameters  for  eq  3. 


Total  water 


content 

(%  of  dry 
weight) 

A 

(%  of  dry 
weight) 

B 

e°C  X  %  of 
dry  weight) 

C 

K 

t°C-cm2/ 

kgf) 

Wu  when  d  *  KS  = 
(%  of  dry  weight) 

to 

1.04 

2.03 

-2.33 

-0.277 

5.77 

12.5 

1.05 

1.94 

-2.45 

-0.292 

5.80 

15 

1.01 

1.96 

-2.34 

-0.236 

5.60 

17.5 

0.905 

2.20 

-2.1  1 

-0.368 

5.55 

20 

0.968 

2.02 

-2.36 

-0.271 

5.74 

22.5 

0.940 

2.32 

-2.18 

-0.339 

6.00 
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Table  2.  Significance  of  surcharge  effect. 


Total  water  content 
(rc  of  dry  weight ) 

Residual 

Surcharge  in  model  mean  square 

1  * 

10 

Yes 

0.0945 

1.41 

1.52 

No 

0.1  333 

1  2.5 

1 .5  2 

1.15 

15 

1.39 

2.00 

17.5 

Yes 

0.0898 

1.68 

2.2  3 

No 

0.1512 

20 

Yes 

0.0748 

1.47 

2.09 

No 

0.1  103 

22.5 

Yes 

0.1098 

1.58 

0.810 

No 

0.1  735 

*  If  F  >  1 .26  the  hypothesis  that  the  variability  or  dispersion  of  the  residuals 
is  greater  for  the  model  when  surcharge  is  not  included  can  be  accepted  at  the 
90%  confidence  level. 

**  If  Z  >  1.28  the  hypothesis  that  the  variability  or  dispersion  of  the  resid¬ 
uals  is  greater  for  the  model  when  surcharge  is  not  included  can  be  accepted 
at  the  90%  confidence  level. 


The  goodness  of  fit  of  eq  3  was  compared  to  that 
of  eq  2  by  analyzing  the  residuals  with  the  nonpar- 
ametric  Siegel-Tukey  test  for  comparing  dispersions 
(Siegel  and  Tukey  1960).  The  statistic  Z,  calculated 
at  each  total  water  content,  is  given  in  Table  2.  For 
four  of  the  values  of  total  water  content,  eq  3  fits 
the  data  better  than  eq  2  at  the  909f  level  of  con¬ 
fidence.  This  claim  cannot  be  made  for  the  other 
two  water  contents. 

We  conclude  that  the  experiments  show  a  signif¬ 
icant  effect  of  surcharge  and  that  the  general  magni¬ 
tude  of  this  effect  is  represented  by  the  values  of  K 
given  in  Table  1 . 

DISCUSSION  AND  CONCLUSIONS 

The  average  value  for  K  in  this  study,  -0.297°C- 
cm2/kg,  is  40  times  greater  than  that  predicted  from 
the  Clausius-Clapeyron  equation  with  surcharge 
pressure  increasing  both  ice  and  liquid  phase  pres¬ 
sures;  it  is  3.4  times  greater  than  predicted  for  sur¬ 
charge  pressure  increasing  only  ice  phase  pressure. 
Williams  (1967)  argued  that  any  surcharge  on  a 
partially  frozen  soil  will  cause  the  pressure  on  both 
the  ice  and  water  phases  to  rise  by  the  same  amount 
at  corresponding  unfrozen  water  content  values. 


This  argument  assumes  that  the  system  has  reached 
thermodynamic  equilibrium.  In  the  experiments 
reported  here  the  soil  samples  were  frozen  first,  and 
then  the  surcharge  was  applied.  Unfrozen  water 
contents  were  measured  at  various  temperatures  as 
the  samples  were  warmed  in  stages  up  to  0°C.  Under 
the  conditions  of  unconfined  loading  used,  it  is  likely 
that  equilibrium  was  not  reached,  but  that  the  sur¬ 
charge  pressure  was  concentrated  in  the  ice  phase, 
causing  the  ice  to  migrate  from  areas  of  high  stress 
concentration  to  areas  of  low  stress  concentration. 
Apparently,  on  the  average,  this  migration  was  not 
completed  during  the  time  the  experiment  was  being 
performed.  If  it  is  assumed  that  the  surcharge  pres¬ 
sure  was  concentrated  by  a  factor  of  3.4  into  the 
ice  phase,  then  the  high  effect  of  surcharge  on  the 
unfrozen  water  content  can  be  explained. 

The  conditions  used  in  this  experiment  correspond 
to  conditions  in  the  field  where  a  surcharge  load  is 
placed  on  an  already  frozen  soil.  The  unfrozen  water 
content  will  be  higher  in  the  stressed  region  than 
what  has  been  predicted  by  the  Clausius-Clapeyron 
equation  using  unfrozen  water  content  data  on  un¬ 
loaded  samples.  This  higher  unfrozen  water  content 
would  cause  the  ice  to  migrate  away  from  the  area 
of  high  stress  faster  than  might  otherwise  be  pre¬ 
dicted. 
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